Fluctuations of Conductance Peak Spacings in the Coulomb Blockade Regime: Role 

of Electron-Electron Interaction. 
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We study influence of electron-electron interaction on 
statistics of Coulomb blockade peak spacings in disordered 
quantum dots. It is shown that the interaction combined with 
fluctuations of eigenfunctions of the Fermi sea, enhances the 
peak spacing fluctuations, in accordance with recent exper- 
iments. In addition, account of the spin degrees of freedom 
leads to a pronounced odd-even structure for weak interaction 
(e 2 /e <C vf)', in the opposite case (e 2 /e > vf) this structure 
is washed out. 

PACS numbers: 73.23.Hk, 73.61.-r, 05.45.+b 



Recent experimental studies of chaotic quantum dots 
in the Coulomb blockade regime showed unusually large 
fluctuations of conductance peak spacings. The r. m. s 
fluctuation of the spacing between consecutive peaks was 
found to exceed the mean level spacing A jj]]^]. These 
fluctuations can not be explained from the point of view 
of the standard "charging energy" model (CEM) ||] 
which assumes that the spacing between the consecutive 
peaks is equal to Sn = e 2 /C + S£_N- Here C is the capac- 
itance of the quantum dot, while 6£n is the distance be- 
tween the last filled single-particle level and the previous 
one. Since the charging energy e 2 /C is a no n- fluctuating 
quantity in CEM, the peak spacing fluctuations are solely 
due to the single-particle energies. Assuming the appli- 
cability of the random matrix theory (RMT) one obtains 
Pj (we consider the case of broken time reversal symme- 
try throughout the paper) 

r. m. s.(Sjv) = A(3tt/8- 1) 1/2 , (1) 

that is less than observed experimentally. 

Below we study the fluctuations of peak spacings in 
the Coulomb blockade regime. First we investigate the 
Coulomb interaction in the finite system, and find cor- 
rections to CEM. Then we demonstrate that these cor- 
rections, being combined with fluctuations of eigenfunc- 
tions, lead to enhancement of spacing fluctuations. This 
becomes important for a sufficiently strong interaction 
and may account for experimentally observable effects. 
In addition, we study the effect of spin degree of freedom 
and find a pronounced odd-even structure of conductance 
peaks, which for strong interactions is washed out. 

Recently Sivan et al Q made an attempt to explain 
the large fluctuations of Sn observed in experiment us- 
ing numerical diagonalization of a Coulomb system with 



N ~ 10 electrons. However, they considered the range 
of very strong interaction e 2 /C > Ep, where Ep is the 
Fermi energy, while in the experiment e 2 /C ~ O.QmeV 
and Ep ~ 12meV. Effect of spin has been also recently 
studied by Prus et al || in terms of a disordered Hubbard 
model. 

We consider a diffusive quantum dot (R ^> I, R and I 
being the size of the dot and the mean free path, respec- 
tively) with a large number of electrons. Most of them are 
compensated by a positive background with the density 
p: in addition, there is a number N of excess (uncom- 
pensated) electrons, so that the total charge of the dot is 
iVe. The Hamiltonian of the system is 



H = 



E 
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Here A labels eigenstates of non-interacting system with 
eigenfunctions tp\(r) and energies e\; the correlations of 
these states have been calculated in Ref. ||. Further- 
more, U(r) = Uo(r) = e 2 /er is the Coulomb interaction, 
e being the dielectric constant. The interaction with pos- 
itive background should be added to Eq. (||). 

Random-phase approximation in restricted ge- 
ometry. As is well known, it is not sufficient to consider 
the Coulomb interaction in the first order of perturba- 
tion theory due to its long-range nature. A common 
way to improve it is the random phase approximation 
(RPA). Whereas it is trivial to solve the RPA equation 
in infinite system, restricted geometry makes the situa- 
tion more complicated (cf 0). In the 3D case the RPA 
equation for the effective potential U is 

U(r, r') = U (r - r') - J dr x dr 2 U a {r - n) 

x K(n,r 2 )U(r 2 ,r'), (4) 
where the polarization operator K has the form 
K (r, r') « v (6(r - r') - V' 1 ) = v ]T <(> a (r)<j> a (r'). 

Here v is the density of states, Q (r) is the eigen- 
function of the Laplace operator with eigenvalue — E a : 



1 



4>q = V~ x l 2 is the so-called zero- mode. In addition, the 
function K contains also a random part. The latter gives 
rise to the fluctuations of the charging energy and is for 
a moment ignored. The bare Coulomb potential Uo can 
be expanded as follows: 

U (r - r') = J2(^ 2 /e)E- 1 ^ a (r)^ a (r') 

+V -V 2 «oa [Mr) + Mr')] + V^uoo, (5) 

where we have denoted 

e 2 f e 2 f 

u ° a = e yi/2 / dr M r )^(r), u 00 = — / dr$(r), 

$(r) = J dr'\r - r'\~ l . (6) 

Then the solution to the equation (0) reads as 

U(r,r') = Y,(^e 2 /e)(E a + n^M^Mr') 

a=£0 

+V~ 1/2 u o a E a (E a + k 2 )- 1 [Mr) + Mr')] 
+V- l (u m -vY,^ + K 2 /E a )- 1 u 2 a ), (7) 

with k = Aire 2 v / e being the inverse screening length. 
One can show Q that in the limit kR ^> 1 (i? is charac- 
teristic size of the system) , which we assume from now on, 
RPA result ([?]) can be also obtained from the Thomas- 
Fermi approximation (TFA), that assumes proportional- 
ity between effective potential and excess charge density. 

Up to a constant, the first term in Eq. (Q) is the usual 
3D screened Coulomb interaction, 

U K {r l7 r 2 ) = (e 2 /er)exp(-Kr), 

where r = |ri — r%\. Two other terms appear due to the 
restricted geometry. The last one is a constant, e 2 /C, 
and C is usual electrostatic capacitance of the system; 
for the sphere with radius R we recover C = eR. Rela- 
tive correction to C due to the finite value of screening 
parameter k is found to be of the order (kR) -1 . Finally, 
the second term (to be denoted as U(r) + U(r')) is a 
single-particle (rather than a two-particle) contribution. 
It vanishes in the electrostatic limit k — > oo, and for the 
sphere takes the explicit form (up to a constant) 

U(r) = -(e 2 /enR 2 ) exp(-n(R - r)), (8) 

i.e. it is localized in a narrow layer of thickness k _1 
near the boundary. Thus, the Coulomb interaction in 
the restricted geometry is given by 

U{r, r') = U K {r, r') + U(r) + U(r') + e 2 /C; 

J drU K (r 7 r') = J drU{r) = 0. (9) 



The form (|Sj) can be easily explained. The static 
screening is the dressing of the added electron by a pos- 
itive cloud. This cloud is created by rearrangement of 
other electrons, and due to the charge conservation an 
excess negative charge e is generated. In the infinite sys- 
tem this excess charge moves to infinity, and does not 
play any role. However, in the finite system it can only 
move to the boundary, and creates the additional con- 
stant potential e 2 /C inside the system. The potential 
U is the edge effect: density of the excess charge grows 
towards the boundary, which in TFA gives rise to the 
additional nonuniform potential U. Thus, the total po- 
tential is the sum of the contributions from the positive 
cloud (U K ) and from the negative charge pushed out to 
the boundary (U + e 2 /C). Note that the terms U K and 
U, ignored by CEM, turn out to be important for the 
peak spacing fluctuations. 

This consideration can be generalized to the 2D geom- 
etry The function K in Eq. ^ has a form (r = (x,y)) 

K(r, r', z, z') = vJ2 Mr)Mr')S(z)S(z'). 

Technically it is convenient to solve this equation in the 
cylinder, and then consider the potential U(r, r') con- 
fined to the plane z = z' = 0. The result has the same 
form (|^) with k = 2ire 2 v/e, the charging energy 

e 2 /C = V-^uoo - v J2(l + KE- l ' 2 )- l ul a \ 

and the 2D screened Coulomb interaction U K (r, r') given 
by 

U K (r,r') = J2(^ 2 /e)(E 1 J 2 + n^M^Mr')- (10) 
The single-particle potential U{r) takes the form 

u(r) = v- 1 ' 2 ]T eH 2 (eH 2 + ^yWoMr) (ii) 

Here Uo a and uqq are given by Eq. (|^) , with the integra- 
tion going over the 2D coordinates. In the limit kR 3> 1 
these results are again equivalent to TFA. In particular, 
for a circle with radius R the potential U (r) is 

U(r) = -e 2 (2nR)- 1 (R 2 -r 2 )- 1 ' 2 . (12) 

Note that, in contrast to the 3D case, U(r) is a smooth 
function of magnitude ~ e 2 / kR 2 all over the sample. 

In the further consideration we replace the bare 
Coulomb potential Uq{t\ — r 2 ) in Eq. (||) by the screened 
potential U{r\,r^) given by Eq. (JsJ) - 

Peak positions. The position of the conductance 
peak fipf in the Coulomb blockade regime is given by 
the variation of the total energy of the system E(N) 
when an additional (excess) electron is added: /ijv = 



2 



E(N + 1) — E(N); the distance between two adjacent 
peaks is the difference Sn = A*jv+i — Mw- 

We use Hartree-Fock (HF) approximation and account 
for e — e interaction by introducing the effective single- 
particle Hamiltonian Hn which depends explicitly on the 
number of excess electrons N. Its spectrum is schemat- 
ically shown on Fig. 1. First empty state is separated 
from the last filled one by a gap of order e 2 /C. The key 
point for the further discussion is that the Hamiltonian 
Hn is essentially random for each N. Namely, we assume 
that the statistical properties of its single-particle excited 
states (with the exception of the gap mentioned above) 
are the same as that of single electron states in a random 
potential. In particular, in the leading approximation 
they obey RMT, an assumption which is in agreement 
with the experiment |10| . 
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FIG. 1. Energy states of the HF Hamiltonian with N (left) 
and N + 1 (right) excess electrons. 

The distance between two adjacent peaks can be de- 
composed in a following way (Fig. 1): 



C _ /..(W+2) (N+2). , (N+2) (JV+lK _ p 
- IMjv+1 ~ J + i^N ~ l*N ) = E\ 



E 2 



Here fi^' is the energy of single-particle eigenstate j of 
the (HF) Hamiltonian Hi. 

Both quantities E\ and E 2 are random. The statistical 
properties of the latter one are trivial, since it is just 
a distance between two adjacent single-particle levels of 
the same random Hamiltonian Hn- Hence, they obey 
RMT; in particular, the average (E2) = A, while the 
fluctuations are given by Eq. (|l|). 

On the other hand, the quantity E\ is the shift of the 
N + 2-th single-particle level due to addition of a new 
electron in the N + 1-th state. It can be expressed as 

Ei = J dr 1 dr 2 U(r ll r 2 ) [\ipN+i(ri)ipN+2{r 2 )\ 2 

- ^N+i{ r i)^N+i{r 2 )^N+2{r 1 )^ N+2 {r2)] . (13) 

Here ipN+i and ipN+2 are eigenfunctions of two first ex- 
cited single-particle states of the Hamiltonian Hn (with 



corresponding energies A t j^ +1 ^ and ^ +2 \ respectively; 
see Fig. 1). The quantities Ei and E 2 are statistically in- 
dependent; therefore fluctuations of Sn are given by the 
sum of fluctuations of E\ and E 2 . In other words, since 
E\ is related to the change in the Hamiltonian (from Hn 
to Hn+i), it leads to fluctuations, additional to those 
given by RMT. 

Statistical properties of this quantity can be investi- 
gated with the use of the expressions for correlations of 
eigenfunctions Q: 

V 2 (\Mri)Mr2)\ 2 ) - 1 = k d (r)[6 k i + n(n,n)] 
+n(n, r 2 )S kl + (i/2)n 2 (n, r 2 ). (14) 

The short-range function kd(r) is defined as 



k d (r) = (nv)- 2 [lmG»]' 



exp(— r/l) 



2D 



(pFr)~ 2 sin 2 pFf, 3D 



and n(ri,r2) is the diffusion propagator. For r 3> I and 
for the case of spherical (3D) and circular (2D) particle 



n(ri,r 2 ) 



{irg)- 1 In R/r, 2D 
2R/(3gr), 3D 



(15) 



with g — 2wE c /A being the dimensionless conductance; 
E c = D/R 2 and D are the Thouless energy and the dif- 
fusion coefficient, respectively. For a moment we assume 
electrons to be spinless, then the wave functions tpN+i 
and tpN+2 are different. 

It can be shown 0], that the average (E{) = e 2 /C, with 
negligible corrections of the order A(ppl) 1 ~ d . Thus, the 
average spacing between the Coulomb blockade peaks is 
(Sn) = e 2 /C + A, in agreement with CEM. 

Now we turn to fluctuations of the quantity E\. In 
accordance with decomposition (||), there are three con- 
tributions. First, the charging energy e 2 /C is a fluc- 
tuating quantity due to fluctuations in the polarization 
operator K(r,r') ||. A direct calculation || shows 
r.m.s (e 2 /C) ~ a 2 Ahig/g. Here the parameters = 
4:Tre 2 v j en 2 and a-i = 2ite 2 v/tK are equal to unity for the 
case of weak interaction (k <C Pf ~ see below). 

Next, we evaluate the fluctuations due to the screened 
interaction U K . Expanding the average of eight wave 
functions in the cumulants jll[] and using Eq.(|l4|) we 
find that these fluctuations are of order r.m.s [Ei(U K )] ~ 
a ( jg~ 1 A. Finally, the fluctuations due to the potential 
U can be directly evaluated with the use of Eq. Jl4|), 
yielding 

r.m.s [Et{U)} 

= [2 ( dr 1 rfr 2 L>(r 1 )L>(r 2 )(|^ 2 r+1 (r 1 )^ 2 r+1 (r 2 )|)] 1 / 2 



(16) 



a 2 g- 1/2 A, 2D, 
«35" 1/2 A, 3D, 
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which constitutes the main contribution to the fluctua- 
tions of the quantity E\ . This result reflects the fact that 
fluctuations of the density of the last added electron lead 
to fluctuations of its energy in the non-uniform potential 
U(r). 

In the above treatment the screened Coulomb poten- 
tial has been calculated within RPA, and thus we have 
assumed k « pj?, that implies ad — 1. If the interac- 
tion is strong (but still the ground state of the dot is 
Fermi- liquid-like), the range of the screened potential is 
determined by the Fermi wavelength. The quantity ad 
is then of order of ad ~ e 2 /evp > 1- The limitation on 
ad given by the Wigner crystallization is a < a c , where 
a c ~ 100 in 2D (see e.g. [jf2)). Therefore, even for a well 
below the Wigner crystallization threshold, the fluctua- 
tions (|l^) can exceed A, if the dimensionless conductance 
g is not too large. 

Although the above results are derived for diffusive 
systems (R ^> I), we believe also that they are valid 
for chaotic ballistic systems (R < I) as well. For the 
latter case the parameter g can be roughly estimated as 
the ratio of inverse time of flight tj to the level spacing, 
g ~ (tfA)^ 1 ~ vp/(RA), up to a geometry-dependent 
coefficient. 

Spin effects. We denote as f and J. two states with 
the same energy but different values of spin; their eigen- 
functions are identical. If the state | is occupied by an 
electron, the energy of the state | is shifted, the shift 
being 



Ei 



drxdrzU^r^xin^xirz)?. (17) 



Using Eq. mj), we find 



(Ei) 



CdAin/pff- 1 Hp f /k), e 2 /e < v F 
~ «dA, e 2 /e>v F 



(18) 



with the coefficients C3 = 1/2 and C2 = 1/ir. In the weak 
interaction regime (e 2 /e <C Vp) we obtain (E\) <C A. 
This means that the set of peaks is split into pairs; the re- 
duced spacing Sn — e 2 /C between two peaks of the same 
pair (states | and |) is equal to (Ei), while that between 
two different pairs is of order A. Furthermore, since the 
eigenfunctions of states f and J, are identical, the corre- 
sponding peak heights are correlated, and the resulting 
picture is a set of pairs with small reduced spacings and 
correlated heights - a pronounced odd-even structure. 
Fluctuations of E\ are given by Eq. (16) multiplied by a 
factor of 2, and do not destroy the pair structure. 

For strong interaction (e 2 /e > vp) the splitting JTs| ) 
exceeds A, and reordering of energy levels takes place. 
However, the spacing fluctuations are still enhanced in 
comparison with the spinless value, since a new type of 
ensemble is created, with the spectrum formed by a su- 
perposition of a RMT-type set of levels and the same set 
shifted by E\ , 

v(E) = ME - Ex) + S(E + Ex- Ex)] . 



In particular, the two-point correlation function is given 
by (A is now the mean level spacing for a system with 
removed spin degeneracy, s — hoj/2A, s — irE±/2A) 

R{oj) = (u(E)u(E + oj))/(v{E))(v{E + u)) 

= ^Rrmt(s) + -Rrmt(s - 5) + jRrmt(s + s), (19) 



where Rrmt(s) = 1 



' sin s is the RMT two-point 



correlation function. The level repulsion in the ensemble 
( |l9| ) is clearly reduced. Thus, fluctuations of peak spac- 
ings are enhanced by spin effects, in addition to their 
enhancement due to fluctuations in E\ discussed above, 
Eq. (|i~6|) . On the other hand, even-odd correlations of 
heights of adjacent peaks in the strong interaction regime 
are washed out. 

In conclusion, we have studied Coulomb interaction ef- 
fects on the statistics of conductance peaks spacings in 
the Coulomb blockade regime. The Coulomb interaction 
leads, in combination with fluctuations of eigenfunctions, 
to enhancement of peak spacing fluctuations in compar- 
ison with the RMT value. Taking into account the spin 
degrees of freedom, we find in addition in the case of weak 
interaction (e 2 /e <C Vp) a pronounced odd-even struc- 
ture. The peaks come in pairs, with correlated heights 
and small reduced spacings S — e 2 /C <C A within each 
pair. In the opposite case e 2 /e > Vp , which is relevant 
to the experiment pL||, this structure is destroyed. 
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